An ultra-broadband electromagnetically indefinite medium formed by aligned carbon nanotubes 
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Anisotropic materials with different signs of components of thie permittivity tensor are called indefinite mate- 
rials. Known realizations of indefinite media suffer of high absorption losses. We show that periodic aiTays of 
parallel carbon nanotubes (CNTs) can behave as a low-loss indefinite medium in the infrared range. We show 
that a finite-thickness slab of CNTs supports the propagation of backward waves with small attenuation in an 
ultra-broad frequency band. In prospective, CNT arrays can be used for subwavelength focusing and detection, 
enhancing the radiation efficiency of small sources. 
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Metal-dielectric nanostructures [1, 2] have already brought 
applications to microscopy and sensing [3] and have the po- 
tential for realization of future nanosized optical devices and 
circuits [4]. In 2003 the authors of [5] noticed the potential 
of so-called indefinite metamaterials for subwavelength imag- 
ing of objects at electrically large distances from them (the 
concept of such imaging was first introduced in [6]). Indef- 
inite metamaterials are artificial uniaxial materials in which 
the axial and tangential components of the permittivity and 
permeability tensors have different signs. In these materials 
the isofrequency surfaces have hyperbolic shape. This re- 
sults in a possibility to design a "hyperlens", where evanescent 
near fields are transformed into propagating modes and can be 
transported at electrically long distances [7]. 

The main challenge on the way to realize this and other ef- 
fects is in the realization of low-loss indefinite materials. For 
waves of only one linear polarization (TM-polarized waves 
with respect to the axis of positive permittivity) it is enough 
to realize a layer of an indefinite dielectric metamaterial [8], 
whose permeabiUty is unity and only components of the per- 
mittivity tensor have different signs. For the visible range such 
a metamaterial was designed in [8] as an array of parallel plas- 
monic (metal) nanowires. Since all plasmonic phenomena are 
related to strong dissipation, the axial component of the per- 
mittivity has a significant imaginary part which strongly re- 
stricts the subwavelength imaging property of the hyperlens 
[8]. In the microwave range, materials with negative permit- 
tivity can be realized as arrays of thin metal wires [9, 10]. 
However, if the field varies along the wires, the properties of 
the structure are more complicated that those of a continuous 
medium due to strong spatial dispersion [11]. The spatial dis- 
persion can be suppressed partially [12] or even totally [13- 
15]. However, manufacturing of structures [12-15] becomes 
a challenge already at millimeter waves as dimensions need 
to be quite small for high-frequency applications. For elec- 
tromagnetic waves in the THz and mid infrared (MIR) range 
there are no known structures with an indefinite permittivity 
tensor. The imaginary and real parts of complex permittivity 
of metals in this range are of the same order and rather high. 

Here we show that arrays of metallic carbon nanotubes 
(CNT) behave as ultra wide-band and low-loss indefinite ma- 
terials in the THz and MIR ranges. This collective property of 
arrays results from the quantum properties of individual CNT, 



namely their high quantum inductance, so-called kinetic in- 
ductance [16], which leads to suppression of spatial disper- 
sion in arrays of parallel CNTs. This can open a possibility 
for realization of interesting devices, which we discuss here. 

Two-dimensional periodic arrays of carbon nanotubes are 
fabricated by many research groups. Such CNT arrays form 
finite-thickness slabs which are used already as field emitters 
[17], biosensors [18] and antennas [19, 20]. We assume that 
all nanotubes possess metallic properties that is a realistic as- 
sumption in view of recent studies of single-wall CNTs [21- 
23]. Usually carbon nanotubes form hexagonal lattices in ar- 
rays in processes of fabrication. However, for low-density ar- 
rays for which arrangement of nanotubes is not important. By 
this reason and for simplicity we assume lattices to be square 
with the constant d. 

For eigenwaves propagating in arrays of infinitely long car- 
bon nanotubes we take the space-time dependence of fields 
and currents as exp [j(wt — k^z)]. Effective electromagnetic 
properties of CNT arrays can be understood from the nonlocal 
quasistatic model of wire media (WM) [15]. In the framework 
of the effective medium model the CNT array can be consid- 
ered as a uniaxial material with the permittivity dyadic 
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where eq is the permittivity of vacuum (we consider CNTs 
placed in vacuum). As was shown in [15], 
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where kp is the effective plasma wavenumber = 
icntC'cnt/(eoMo)^ Lent, Cent are the effective inductance and 
capacitance per unit length, respectively, and the parameter 
^ = (i?cnt/icnt)v^eoMo is responsible for losses. Distributed 
parameters Lent, Cent and _Rcnt for a separate carbon nanotube 
can be obtained using the model of impedance cylinder and ef- 
fective boundary conditions developed in [24] . According to 
this model the simple approximate expression for the complex 
surface conductivity, which is valid for metallic zigzag CNTs 
at frequencies below optical transitions looks as: 
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Here e is the electron charge, Fq = 2.7 eV is the overlap- 
ping integral, t = 1/j^ is the relaxation time, h is the reduced 
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Planck constant. The radius of metallic zigzag CNT, taken as 
an example, r is determined by an integer qasr ~ 3V3qb/2n, 
where b = 0.142nm is the interatomic distance in graphene. 
Since the wall of a single-wall CNT is a monoatomic sheet of 
carbon, (3) can be considered as the surface conductivity of 
the carbon nanotube. The surface impedance per unit length 
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= RQ+jujLo. (4) 



Let us find Lent, Cent and i?cnt for a separate carbon nan- 
otube. The total inductance per unit length Lent — Lq + Lem, 
where Lo is the kinetic inductance [16] defined by the formula 
(4), which has a quantum nature, and Lcm is the electromag- 
netic inductance per unit length. For thin CNTs the kinetic in- 
ductance strongly dominates over the electromagnetic induc- 
tance, e.g., for the zigzag CNT, having the radius r ~ 1.53 nm 
(q ^ 13), Lom = 5.6 X lO"'^ H/m and Lq = 3.7 x lO^^ H/m. 
This shows that the electromagnetic inductance can be ne- 
glected. The total capacitance of CNTs Cent includes the elec- 
trostatic capacitance defined as Cs = eo/ito/Lem and the quan- 
tum capacitance Cg = 2e^ /{hvp), where vp is the Fermi 
velocity which is equal to 8 x 10^ F/m for graphene and 
CNT [25]. The two effective capacitances are connected se- 
rially. For the considered example Cs — 1.98 x 10^^^ F/m, 
Cq ~ 9.66 X 10-" F/m and Ctot = CeCg/{Ce + Cg) = 
1.64 X lO^^^F/m. Let us estimate the terms entering for- 
mula (2) at the frequency uj /2Tr = ITHz. One obtains 
= 4.32 X IQi" l/m2, P = 4.39 x 10*^ l/m^ and dimen- 
sionless — 4.32 x 10^, so the last term in the denominator 
of formula (2) can be neglected. 

Thus strong spatial dispersion is suppressed and the mate- 
rial behaves as a uniaxial free-electron plasma, where elec- 
trons can move only along z-direction. Apparently, this cor- 
responds to an indefinite medium at frequencies below the 
plasma frequency, because the permittivity in the directions 
orthogonal to the tubes is close to unity. Dispersion equation 
for waves propagating in a uniaxial crystal is [26] 



where k'j_ = k^ + ky. Substituting (2) with fc^/n^ 
Eq. (5) we obtain a simple formula: 



k\e 



kl 



k'^-kl 



(5) 



into 
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which describes a typical conic-type dispersion, see Fig. 1. 
Applicability of the described effective medium model was 
verified comparing with the results of the electrodynamic 
model, based on Green's function method [27]. This model 
takes into account the array periodicity and demonstrates that 
the effective medium model gives a good agreement with 
the more accurate theory for small transversal wave numbers 
{\k±d\ < O.Wtt). 

In a series of papers published in 2003-2004 by the groups 
of G. Eleftheriades, C. Caloz and T. Itoh, Pendry's concept of 
perfect lens [6] was experimentally confirmed for a planar ana- 
logue of the perfect lens (see also books [28] and [29]). Pla- 
nar networks composed of backward-wave transmission lines 




FIG. 1 : (color online). The real part of the slow-wave factor in the 
transversal plane, calculated at ijj/2n = 1 THz for d — 2.5 nm. 




FIG. 2: (color online). Frequency as a function of the normalized 
wavenumbers in the transverse plane kxd/n, kyd/n. d — 10 nm, 
and h — 1.5 fim. 



(TLN) are surface analogues of bulk double-negative mate- 
rials. Surface formed by TLN acts as a layer of an indefi- 
nite medium, where the spatial dispersion is suppressed [30]. 
Based on the results of the previous duscussion, we expect that 
in the THz and MIR ranges a similar backward-wave structure 
can be realized as an array of aligned CNTs. 

Let us consider waves propagating in a finite-thickness slab 
of CNTs assuming for simplicity that the slab is placed be- 
tween perfect electric conductor (PEC) and perfect magnetic 
conductor (PMC) planes, where the PMC boundary models 
the open-ended interface with free space. The resonant condi- 
tion for the longitudinal wavenumber k^ reads kz = 7r/(2/i) 
where h is the thickness of the slab. The relation between the 
transversal wave vector components and the wavenumber in 
free space is the following: 



k I — k^ 



ky - 



k^ 



(7) 



One can easily show that the derivative dfc^/dfc^ < if 
kz/k > 1 and kp/k > 1. Thus, in this regime the finite- 
thickness slab supports propagation of backward waves, be- 
cause the group velocity is in the opposite direction to the 
phase velocity. 

This property can be understood also from considering a 
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planar waveguide filled with a CNT array (the array axis is 
orthogonal to the walls of the waveguide). The propagation 
constant along the waveguide is equal to 

= Ve.. - (m7r/2/i)2], (8) 

where m is a positive integer and h is the height of the waveg- 
uide [11]. If ezz < 0, backward- wave propagation is allowed 
when k < rmr/2h and forbidden for k > nnr/2h. 

Fig. 2 illustrates dispersion properties of three modes, corre- 
sponding to m = 1,2 and rn = 3. There are three embedded 
conical surfaces. The internal cone corresponds to to = 1, 
and the external one to m = 3. One can see that backward 
waves propagate in the slab in a very wide frequency range. 
Their properties are quite isotropic in the xy-plane due to a 
very small period of the CNT lattice. It follows from formula 
(8) that there is no low-frequency cutoff for any mode, which 
can propagate at low frequencies with very large transversal 
wavenumbers k± . The effective medium model becomes un- 
applicable if k±d > O.IStt (compare with results of [27]). 

To suppress the spatial dispersion we do not need any 
loadings because CNTs possess very high kinetic inductance 
which ensures that k1/n^ <C k^- Comparing to the microwave 
backward-wave structures based on loaded transmission lines 
and mushroom layers we note that in CNT arrays there is no 
parasitic inductance between the cells and parasitic capaci- 
tance between the tube ends and the ground plane is very small 
due to the small cross-section area of the tubes. This leads to 
a dramatically increased bandwidth of backward-wave propa- 
gation and gives the ground to consider CNT arrays as perfect 
backward-wave metamaterials. 

For any optical material at terahertz and infrared ranges 
one of the most important issues is the level of losses. In 
CNTs it is determined by the relaxation time t at the fre- 
quencies below optical transitions. At low frequencies, in- 
cluding values close to a few terahertz, t at room temperature 
usually is estimated as 3 x 10'^^ s [24, 31]. Such a relax- 
ation time gives the ratio Im(fcj^)/Re(A:^) ~ 10^'^ at frequen- 
cies around 50 THz. However, in other sources this value is 
taken to be 10^^'^ s [32], which corresponds to considerably 
higher losses, namely, Im(fcj^)/Re(fcj^) ~ 0.025 in the same 
frequency range. 

Here we point out one potential application of indefinite di- 
electric materials based on the conversion of incident evanes- 
cent waves into propagating transmitted waves by samples of 
an indefinite medium. Such a conversion is possible up to very 
high, in the limiting case infinite spatial frequencies. The idea 
is illustrated by Fig. 3. 

Consider a prism with the wedge angle 2a « tt/2 filled 
with an indefinite medium with components of the permittiv- 
ity tensor e^z < 0, £yy = £xx > 0. For a wedge made of a 
CNT array this would correspond to the z-oriented tubes. Let 
an evanescent wave attenuating along z and having transversal 
wave vector component ky > k = LUy/eoJIo impinges on the 
left side of the prism. In the isofrequency diagram {k^ — ky) 
which is depicted in Fig. 3 (a) the wave vector of the refracted 
wave selects the point at the medium isofrequency contour so 
that the component of the refracted wave vector which is par- 
allel to the first asymptote of the hyperbolic isofrequency con- 
tour preserves. On the second side of the prism the compo- 
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FIG. 3: (color online), (a) - Conversion of evanescent waves with 
arbitrary spatial frequencies into propagating ones by an indefinite 
wedge with the wedge angle 2a ~ n/2 and horizontal optical axis, 
(b) - Enhanced emission of a poor emitter from a cut in a shell of 
indefinite medium with radial directions of the local optical axes. 



nent of the refracted wave vector which is parallel to the sec- 
ond asymptote of the hyperbola preserves. This holds because 
both sides are tilted with respect to axis y under the same angle 
a as the asymptotes of the hyperbolic isofrequency contour. 
One can easily derive for this geometry that all incident waves 
with spatial frequencies ky < k/ \ cosacos2Q;| transmit into 
propagating ones (ky < k) behind its second interface. If 
a = 7r/4, all evanescent waves are transformed into propagat- 
ing ones. The isofrequency contour with a = tt/A holds for 
medium of parallel CNT at fc = kp/ ^/2. The peculiarity of the 
suggested structure compared to the well known Kretschmann 
dielectric prism widely used in optics for coupling evanescent 
optical devices (waveguides and resonators) to free space is 
the possibility to convert evanescent waves with arbitrary spa- 
tial frequencies at a certain frequency. In the Kretschmann 
method the maximal spatial frequency where the conversion 
of the evanescent wave into propagating wave (or reverse) is 
restricted by the permittivity of the material. 

The suggested effect can be practically used for the en- 
hancement of emission from poor emitters, i.e. light sources 
which low radiation but high level of the electromagnetic en- 
ergy stored in them. The poor emitter can be placed in a spher- 
ical shell of an indefinite material with space-variable optical 
axis. Namely, the local optical axis is directed radially (it can 
be implemented as an array of radially oriented CNT). If the 
shell is spherical, the evanescent waves convert on its internal 
surface to propagating waves and totally internally reflect on 
the outer interface. Now we make a cut in the shell as it is 
shown in Fig. 3 (b). In this part of the shell the inner and outer 
boundaries of the indefinite medium are strongly not parallel 
and one can obtain the radiation dramatically enhanced com- 
pared to the original radiation of the emitter 

We have shown that arrays of single-wall metallic carbon 
nanotubes behave as indefinite media. Such properties are 
provided by a very high kinetic inductance of thin carbon nan- 
otubes. Arrays of finite-length CNTs support propagation of 
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backward waves, which are characterized by low levels of 
losses in the terahertz and mid-infrared ranges. The finite- 
thickness carbon nanotube array can be considered as the per- 
fect isotropic backward-wave metamaterial. We have theoret- 
ically demonstrated that these materials can be used as trans- 
formers of full spectrum of evanescent electromagnetic waves 
into propagating ones, and as devices for emission enhance- 
ment for sub wavelength- sized radiators. 

This work has been partially funded by the Academy of Fin- 
land and Nokia through the Center-of-Excellence program. 
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